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Hamiltonian systems can be classified into ten classes, in terms of the presence or absence of time- 
reversal symmetry, particle- hole symmetry and sublattice/chiral symmetry. We construct a quantum 
coherent scattering theory of linear transport for coupled electric, heat and spin transport; including 
the effect of Andreev reflection from superconductors. We derive a complete list of the Onsager 
reciprocity relations between transport coefficients for coupled electric, spin, thermoelectric and spin 
caloritronic effects. We apply these to all ten symmetry classes, paying special attention to specific 
additional relations that follow from the combination of symmetries, beyond microreversibility. We 
discuss these relations in several illustrative situations. We show the reciprocity between spin-Hall 
and inverse spin-Hall effects, and the reciprocity between spin-injection and magnetoelectric spin 
currents. We discuss the symmetry and reciprocity relations of Seebeck, Peltier, spin-Seebeck and 
spin-Peltier effects in systems with and without coupling to superconductors. 
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I. INTRODUCTION 

Onsager's reciprocity relations are cornerstones of 
nonequilibrium statistical niechanicJ32]_ They relate lin- 
ear response coefficients between flux densities and ther- 
modynamic forces to one another. They are based on 
the fundamental principle of microreversibility which, 
for systems with time- reversal symmetry (TRS), says 
that "if the velocities of all the particles present are 
reversed simultaneously the particles will retrace their 
former paths, reversing the entire succession of configu- 
rations'Cl. When TRS is broken, microreversibility fur- 
ther requires to invert all TRS breaking fields, which, 
to fix ideas, one may take as magnetic fields, fiuxes or 
exchange fields. Combining all of them into a single 
multi-cor npo nent field 5£, the Onsager reciprocity rela- 
tions reacP^ 

£,,(?£) =£,,(-5{), (1) 

where the linear coefficient Lij determines the response 
of the flux density 3i - for instance the electric or heat 
current - to a weak thermodynamic force Xj - for in- 
stance an electric field or a temperature gradient. Thus 
the precise form of the Onsager reciprocity relations de- 
pends on the symmetries of the system. Seminal works 
have classified noninteracting quantum mechanical sys- 
tems into ten general symmetry classe^SHU^ and it is the 
purpose of the present manuscript to derive Onsager's 
relations for all these symmetry classes. Four of them, 
in particular, combine two different types of quasipar- 
ticleJ^, with microscopic representations including, e.g. 
hybrid systems where quantum coherent normal metal- 
lic conductors are connected to superconductors. Since 
at sub-gap energies an interface between a normal metal 
and a superconductor blocks heat currents but not elec- 



tric currents it is natural to ask whether the Onsager 
reciprocity relation between, say, the Seebeck and Peltier 
thermoelectric coefficients survives in such systems. On- 
sager relations in the presence of superconductivity have 
been discussed rather incompletely until novP^^, despite 
much experimental"^**-'*' and theoreticaPtSIIIlill interest 
in thermoelectric transport properties of hybrid normal- 
metallic /superconducting systems. 

Further motivation is provided by fundamental aspects 
of spintronic^^ and spin caloritronics (spin-Seebeck and 
spin-Peltier effects^!, where Onsager relations are of sig- 
nificant interesP^'^Sl. As a matter of fact, reciprocity 
relations decisively helped in experimentally uncovering 
elusive spin effects, by suggesting to measure electric ef- 
fects that are reciprocal to them. As but one example, we 
mention the inverse spin Hall effectPsJ^, where a trans- 
verse electric current or voltage is generated by an in- 
jected spin currenlpS. Onsager relations als o put con- 
straints on the measurement of spin currentJ^^Ell that 
can be circumvented in the nonlinear regime 
Accordingly, we will incorporate spin currents and ac- 
cumulations into our formalism. Several of the Onsager 
relations for spin transport we derive below appeared in 
one way or another in earlier publications, see in par- 
ticular Refs. [2Tll26|32] . Here, we summarize them in a 
unified way and extend them to all ten symmetry classes. 
We are unaware of earlier discussions of Onsager relations 
for spin transport in the presence of superconductivity. 

The classification into ten different symmetry classes 
has recently received renewed attention, because the ex- 
istence of topologically nontrivial phase^^H depends on 
the system's symmetries and its dimensionalitjl^?'. The 
Onsager relations we derive below depend only on funda- 
mental symmetries and are eq ually valid in topologically 
trivial and nontrivial statea^Sl. In several instances, how- 
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Symmetry class 


TRS 


PHS 


SLS 


Physical example 


Wigner-Dyson 


A (unitary) 











mag. flux (no SC) 




AI (orthog.) 


+1 








no mag. flux & no spin-orbit (no SC) 




All (sympl.) 


-1 








spin-orbit & no mag. flux (no SC) 


Chiral 


AIII (unitary) 








1 


mag. flux & bipartite lattice (no SC) 




DUL t^ortnog.j 


-1-1 


-1-1 


1 


nn Tnflo' flux Tin '^niTi-nrVtit f^iT^flrtitp Ifittipp Inn Sir^i 




CII (sympl.) 


— 1 


-1 


1 


spin-orbit & no mag. flux & bipartite lattice (no SC) 


Altland-Zirnbauer 


D 





+1 





SC, mag. flux, & spin-orbit 




C 





-1 





SC, mag. flux, & no spin-orbit 




DIII 


-1 


+1 


1 


SC, no mag. flux, & spin-orbit 




CI 


+1 


-1 


1 


SC, no mag. flux, & no spin-orbit 



Table I: The ten-fold symmetry classification of Hamiltonians. The second column from the left refers to Cartan's nomenclature 
for symmetric spacei^, the three middle columns indicate whether the classes have broken (0) or unbroken (±1) time-reversal 
symmetry (TRS), particle-hole symmetry (PHS) and sublattice symmetry (SLS), while the rightmost column mentions mi- 
croscopic realizations in each class, with SC indicating the presence of superconductivity. Aside from the indicated bipartite 
lattice models with preserved sublattice symmetry, the chiral classes are also realized in low-energy models for quantum chro- 
modynamicp'. Spin rotational symmetry (SRS) is present in classes AI, BDI, C and CI, absent in classes All, CII, D and DIII, 
and irrelevant in classes A and AIII. 



ever, specific additional relations exist, that arise because 
of the conservation of each quasiparticle species (relevant 
to systems without superconductor so that there are no 
Andreev processes converting electrons into holes and 
vice- versa), the presence of particle- hole symmetry or 
sublattice/chiral symmetry. This is the case, for instance, 
for two-terminal thermoelectric transport in the Wigner- 
Dyson symmetry classes where the relation between 
Seebeck, B and Peltier, F coefRcients reads equivalently 
B{K)To = r(J{) or B{^)To = r(-Jf), with the base 
temperature Tq, because of the additional reciprocity re- 
lation r(J£) = T{~3-Cj^. Below, we pay special attention 
to these nongeneric relations. 



The manuscript is organized as follows. In Sec- 
tion [nj we discuss the ten symmetry classes, and the 
crossover between Altland-Zirnbauer and the Wigner- 
Dyson classes as the temperature is raised in hybrid sys- 
tems. In Section |III| we derive and list the symmetries 
that the S-matrix satisfies in all classes. Onsager rela- 
tions will follow from these symmetries, once they are 
inserted into scattering theory expressions for the linear 
transport coefficients. In Section |IV[ we formulate the 
problem in terms of the scattering matrix of the system 
and connect the Onsager coefficients to the system's scat- 
tering matrix. In Section|V]we list the general reciprocity 
relations and mention additional ones occurring in spe- 
cial circumstances. Finally, in Sections |VI| and |VII[ we 
discuss some cases of importance and the associated On- 
sager relations in systems with coupled electric and spin 
transport, superconductors or chiral symmetries. Con- 



II. THE TEN-FOLD WAY 

Hamiltonian systems are classified according to the 
presence or absence of fund ame ntal symmetries. The his- 
torical classification schem^^ is based on TRS and spin- 
rotational symmetry (SRS). Three Wigner-Dyson classes 
are defined in this way. Using the Cartan nomenclature 
for symmetric spaces, the class A has both symmetries 
broken, the class AI has both symmetries present, and 
the class All has broken SRS but unbroken TRS. When 
TRS is broken, the presence or absence of SRS only af- 
fects the size of the Hamiltonian matrix — and not its 
symmetry — and there is thus no fourth class. 

Chiral classes were next introducecP, which cap- 
ture the structure of the QCD Dirac operators. Be- 
side relativistic fermions, they are also appropriate 
to describe bipartite lattice Hamiltonians with unbro- 
ken sublattice symmetry (SLS). Examples include two- 
dimensional square and hexagonal lattices, as well as 
three-dimensional cubic lattices without mass/on-site 
term, the latter generically breaking SLS. Here also, there 
are three classes, with (apart from their chiral symmetry) 
the same symmetries as the Wigner-Dyson classes. 

Finally, four more classes of Bogoliubov-de Gennes 
(BdG) Hamiltonians appear - the Altland-Zirnbauer 
classeJ^ - when normal metals are brought into contact 
with superconductors: with SRS (C and CI) and without 
SRS (D and DIII), with TRS (CI and DIII) and without 
TRS (C and D). When dealing with such systems, we use 
a convention where the BdG Hamiltonian reads 



H 



h~ Us 



A 



(2) 



elusions are given in Section VIII 



A* Use - (T^y'>h*(7(-y'> I ' 

with the Pauli matrix (t'^^ acting on the spin degree 



3 



Symmetry class 


TRS 


PHS 


SLS 


Physical example 


Crossovers 


D ^ A 





+1 ^ 








as D but (rAndrfcaro) is not small 


for Andreev 


C ^ A 





-1 








as C but(rAndrfcB7o) is not small 


interfere. 


Dili ^ All 


-1 


+1 ^ 





1 ^ 


as Dili but (rAndr^BTo) is not small 




CI ^ AI 


+1 


-1 





1^ 


as CI but (rAndrfcBTb) is not small 



Table II: Crossovers from the Altland-Zirnbauer to the Wigner-Dyson classes. Increasing the temperature breaks PHS so that 
the quasiparticle excitation energy e cannot be treated as small. One way to break PHS is to make (fcB7o)''"Andr not negligible, 
where TAndr is a timescale associated with impinging on or returning to the superconducting contacts. 



of freedom and the chemical potential /isc on the su- 
perconductor. With this convention, used for exam- 
ple in Refs. |39|40j . the second-quantized Bogoliubov-de 
Gennes Hamiltonian is ^c^iJc, where 

with being the vector of creation operator for all k- 
states of spin-f electrons, etc. This has the hole sec- 
tor rotated by zcr*^*'^ with respect to the Hamiltonian in 
Refs. |5|37j . The form of Eq. ^ has the advantage that 
upon assuming SRS (so h* commutes with a^^"^), it im- 
mediately reduces to that used in Refs. |7l41H43j . 

Ref. [37] introduced a unifying ten-fold classification 
scheme for all the above Hamiltonians. They considered 
TRS and particle-hole symmetry (PHS), which can both 
be represented by antiunitary operators, and accordingly, 
these two symmetries can be either broken, or unbroken. 
In the former case, we represent this by a 0, while in the 
latter case, the antiunitary operator squares to either 
or —1. A squared TRS of -1-1 corresponds to spinless or 
integer-spin particles, while a squared TRS of —1 corre- 
sponds to half-integer-spin particles. A squared PHS of 
-f 1 corresponds to triplet pairing, while a squared PHS 
of —1 corresponds to singlet pairing in a Bogoliubov-de 
Gennes HamiltoniaiP^. Naively one would think that 
this leads to 3 x 3 = 9 classes, however there are two dis- 
tinct possibilities when both TRS and PHS are broken. 
In this case, the symmetry represented by the product 
of the two antiunitary operators gives either (when the 
corresponding symmetry is broken) or -1-1. This finally 
gives 3x3— 1 + 2 = 10 symmetry classes again. Using 
the just defined three indices, we summarize the ten sym- 
metry classes in Table [l] where we additionally mention 
relevant physical realizations for each of them. 

In this classification, the possible symmetries that the 
Hamiltonian H satisfies are (i) TRS : H = THT-^, 
with T = —iK, with the complex conjugation opera- 
tor K in the spinless case, and T = —ia^^^K for spin-1/2 
fermions, with the Pauli matrix a^^^ acting in spin space; 
(ii) PHS : H = -PHP-\ with P = -ia^y^^v^K, with 
the Pauli matrix r^^^ acting in Nambu spac^^ (iii) SLS 
: H = —ri^^'>Hri^^\ with the Pauli matrix 77'^) acting on 
sublattice space (bipartite lattices are assumed here). 

TRS, SRS and SLS can be broken by an orbital mag- 
netic field, spin-orbit interaction and mass/on-site terms 
respectively. The Altland-Zirnbauer classes assume PHS, 



which strictly speaking forces thermoelectric effects to 
vanish identically. PHS can be broken, for instance, by 
moving away in energy from the special e = symmetry 
point - the superconductor's chemical potential. This oc- 
curs upon increasing the temperature, when the latter ex- 
ceeds a Thouless energy scale, Et ~ ^And,. where the time 
scale TAndr IS related to the time it takes to impinge on 
or return to the normal- metal/superconductor interface. 
This energy is implicitly assumed to be much smaller 
than the superconductor's critical temperature. Table [TT| 
summarizes the crossovers from the Altland-Zirnbauer 
to the Wigner-Dyson classes as the temperature is raised 
such that the coherence between electron and Andreev- 
refiected hole quasiparticle gets lost. We will get back to 
this point in Section [VII C| below. The existence of SLS 
also requires that the spectrum is symmetric about zero 
energy, thus, at half-filling, SLS also leads to the vanish- 
ing of thermoelectric effects, which one recovers as the 
electrochemical potential is tuned away from half-filling. 



III. SYMMETRIES AND RECIPROCITIES OF 
THE S-MATRIX 

Our investigations are based on the scattering the- 
ory of quantum transportP^'^, which, for noninteract- 
ing systems, allows to straightforwardly derive Onsager 
reciprocity relations solely from the symmetries of the 
system's scattering matrix S. 

Reciprocity relations for § follow directly from microre- 
versibilitjB^. They reacP^, 

S(J{) =a(*')S^(-J{)a(^\ (4) 

where a'^y^ is a Pauli matrix acting in spin space, and ' " 
indicates the matrix transpose of spin, transport channel 
and (with superconductivity) quasiparticle indices. In- 
cluded in Eq. Q is the relation S(J{) = §'^(-J{) valid 
when the antiunitary TRS operator squares to 1 and SRS 
is not broken. Eq. Q can be derived by constructing § 
first with scattering states 0„o.(iK), then with their time- 
reversed —ia^y^Kcj^n^^—J-C), with the complex conjuga- 
tion operator K, and equating the two resultJ^. Eq. Q 
is intimately related to Kramers degeneracy, which in 
the presence of TRS (5{ = 0) follows from the symmetry 
property H — a'^y^ H* a'^y'^ of the Hamiltonian H. Speci- 
fying to half-integer-spin particles, it can equivalently be 
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Figure 1: Microreversibility operating on a two-terminal S-matrix in the absence (top) and in tiie presence of superconductivity 
(bottom). The scattering amplitudes a and b depend on spin indices, a =t (+1), 4- (~l)i quasiparticle indices, /i = e(-|-l), 
/i( — 1) and terminal indices 1 and 2. Time-reversal implies inverting the particle's momentum, spin and quasiparticle isospin, 
as well as magnetic fields and fluxes. 



rewritten in a form that renders its connection to mi- 
croreversibility more evident 

where i,j are transport channel indices, ^,,v = e,h are 
quasiparticle indices and a = —a are spin indices. 

Further relations can be constructed by combining 
Eqs. Q and ([s]) with additional symmetries of the §- 
matrix. The latter are obtained by translating PHS and 
SLS of the Hamiltonian into symmetries of the §-matrix. 
For this purpose, we use the relatiorPSl 

B{e) = l + 2TriW\H - e~iTTWW'')-'^W , (6) 

between § and H, with a rectangular matrix W that 
couples the scatterer to external leads. We consider first 
PHS. The presence of superconductivity requires to in- 
troduce electron and hole quasiparticles, and when PHS 
is present, the energy spectrum is symmetric about zero 
energy (taken as the chemical potential of the supercon- 
ductor). With the convention of Eq. ([2]), PHS reads 
H = -a'-y^T'^y^H*a'^v^T'^y'^, with the minus-sign indi- 
cating how the symmetry of the energy spectrum differs 
from Kramers degeneracy. From this we 

obtairP2l 

§(M) = a^y^T^y') §* (j£) a^y)T^y) . (7) 

Combining Eqs. ([4| and ([7| and paying attention to the 
ordering of spin-indices given in Eq. (§, one obtains. 



where quasiparticle indices = +l(e),— l(/i) when 
they appear as pref actors. 

The reciprocity relations ^ and ([7| for the S-matrix 
are illustrated in Fig. [ij We defined the S-matrix via 
the relation § a = b between vectors a and b of compo- 
nents for incoming and outgoing quasiparticle flux am- 
plitudes, respectively. Each component of these vec- 
tors corresponds to a given terminal, a transverse trans- 
port channel in that terminal, a spin orientation and, in 
the presence of superconductivity, a quasiparticle index. 
Complex conjugation in Fig.[I]and Eq.([4]) occurs because 
TRS and PHS are represented by antiunitary operators, 
i.e. products of a unitary operator with complex conju- 
gation. 

We finally comment on SLS. The chiral Hamiltonian 
symmetry reads H = —r]^^^ H r]^^\ with the Pauli matrix 
7](^) acting in sublattice space. For the scattering matrix, 
this translates into 

S(Jf,6)=r?(^)St(J{,-e),7(^), (9) 

where in contrast to earlier symmetry relations, we ex- 
plicitly had to write the energy-dependence of the §- 
matrix. Combining Eqs. ^ and (|9]), one obtains, 

S™",,,(:K,e) = mn(S™?,,(J{,-6))* 
= aa'§p,,^(-5£,e) 
= mnaa'(§^X-5{,-e))*, (10) 

where we introduced sublattice indices m,n = 
A(+1),S(-1). 
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IV. SCATTERING APPROACH TO TRANSPORT AND FORMULATION OF THE PROBLEM 



We consider a multiterminal device connected to i, j = 1, 2, . . . iV electrodes. The linear response relation is 
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(11) 



between heat, J^, electric, /l'^' and spin, /I"-* currents on the one hand, and temperatures, Tj, voltages, Vj and spin 
accumulations, /i^"^ on the other. 



The coefficients with superindices are the usual 
thermoelectric coefficients, while the coefficients G*^"^-* 
are conductances and spin-dependent conductances, re- 
lating electric and spin currents to electric voltages and 
spin accumulations'^"^ . Finally, one has spin-Peltier ma- 
trix elements F^"'^' connecting heat currents to spin accu- 
mulations and spin-Seebeck matrix elements B^"'^^ con- 
necting spin currents to temperature differenced^. The 
resulting spin caloritronic (spin-Seebeck and spin- Peltier) 
effects have been investigated theoreticallySMSSl and ex- 
perimentallji^nmi!^ As usual, we assume that there is no 
spin relaxation in the terminals where spin currents are 
measured, so that the latter are well defined. Our goal is 
to determine reciprocity relations between the elements 
of the Onsager matrix defined on the right-hand side of 
Eq. ^T|) in the ten symmetry classes discussed in Sec- 
tion |nf5 

We will express the matrix elements of the Onsager 
matrix in Eq. ( 11 ) in terms of the S-matrix. We discuss 



separately purely metallic systems and hybrid systems 
consisting of normal metallic components connected to 
superconductors. 



A. Purely metallic systems 

Purely metallic systems fall in either one of the 
Wigner-Dyson or in one of the chiral classes. We start 
from the expression for electric current in Ref. |47|, ex- 
tending it to account for heat and spin currents, e.g. along 
the lines of Refs. [44|53j . This gives us the following lin- 
ear relations between electric, heat and spin currents, on 
one hand, and voltages, temperatures and spin accumu- 
lations, on the other hand; 



de 



dl 
de 



x[Aif)+<5o/3e(T,--To)/To], (12a) 



r(") 



x[^if +SopeiTj-To)/To], (12b) 



where the sums run over all terminal indices z,j and all 
charge-spin indices a, /? = 0, x, y, z. The electrochemical 
potential in terminal j is = fip + eVj with the applied 
voltage Vj and Tq is the base temperature about which 
the Fermi function / = (exp[e/T]-|-l)~-'^ is expanded. The 

spin accumulations \ /? 7^ 0, are one half times the 
/3-components of the spin accumulation vector fij , giving 
the difference in chemical potential between the two spin 
species along the /3 axis, e.g. ^'f^ = {^^p — p!'^'') /2. They 
are nonequilibrium spin accumulations whose origin is of 
little importance here. 



In Eqs. (12), we introduced the spin-dependent trans- 



mission and reflection coefficients 



7)' 



Ir[(S,,)tal"^S,, 



(/3)i 



(13) 



where = 0, x, y, z are Pauli matrices (cr^^-' is the 

identity matrix) and the trace is taken over both spin and 
transmission channel indices. Note the position of the 
Pauli matrices, where cr|"'' measures the spin in direction 
a as the electron exits the systems, while u|^^ m easures it 
along /3 as the electron enters the systenp2IMl_ These co- 
efficients depend on the energy e of the injected electrons, 
which we explicitly wrote in Eq. (12). Reciprocity rela- 



tions for the Onsager matrix elements in purely metallic 
systems directly follow from combining Eqs. ( 12 ) with the 

^"'^^ under microreversibil- 



transformation rules for the T', 
ity. Pauli matrices satisfy ct';^'^^ = (—1) 
nx.y,z — 1 and no = 0. Using this and Eq. ([5| we obtain 

'^"'^^(J{,e) = (-l)"°+"^*Tjf"^(-J£,e) 



(14) 



Thus the reciprocity relation between spin-dependent 
transmission coefficients in Eq. ( 14 ) picks up a minus 



sign if the spin is resolved upon entering the system, and 
another if it is resolved upon leaving the system. 
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B. Metallic systems with chiral symmetry 

The chiral classes correspond to systems with a bipar- 
tite lattice, however currently no experiments are capa- 
ble of measuring sublattice-resolved currents. Thus the 
charge and spin-transport is given by Eqs. ( 12|13 ), with 
the trace over channels supplemented by a trace over the 
sublattice indices (A and B sites). If one could mea- 
sure sublattice isospin current, then one would have to 
add further Pauli matrices acting in sublattice space into 
Eqs. (13 1, leading to extra factors (— l)"o'+"/3' due to 
isospin in Eq. (14). We do not consider this possibility 



further, due to its lack of physical implementation. 

A relevant consequence of SLS is however that from 
Eq. (flo|, we get 7'^f \e) = 7^f^"\-e). Combining this 
with Eq. ( 14 ) we obtain 



= (-1)" 
= (-1)" 



•'ij 



(-3^,-6). (15) 



These relations are strictly valid only insofar as leads are 
preserving SLS, meaning that they connect equally to 
both sublattice sites of each unit cell. 



C. Hybrid superconducting- normal metallic 

systems 

Hybrid normal-metallic/superconducting systems have 
Andreev electron-hole scattering. This scattering may 
induce PHS, in which case the system falls in one of the 
four Altland-Zirnbauer symmetry classes in Table [ip. 

To include Andreev scattering, one has to consider two 
kinds of quasiparticles (electrons and holes), which carry 
excitation energy ±e counted from the chemical poten- 
tial of the superconductor ^sc- These quasiparticles are 
converted into one another when they hit the supercon- 
ductor. Ref. [7] constructed a scattering theory of ther- 
moelectric transport which include these effects. We need 



to include spin currents and accumulations. 

To do this, we go back to the derivation of the scat- 
tering theory in terms of creation and annihilation op- 
erators acting on scattering states in the lead (see e.g. 
Ref. [H]). We write hole creation operators at energy e, 
in terms of electron annihilation operators at energy —e 
S'"/°"(-e) with 



as ci;:)'"/-*^^) 



^(h)in/outt 



(h)in/outt ' 

i;»T 

(h)in/out'i* 



(c)in/out 



^(c)in/out ~ 
(c) in/out 



Here, i gives the index of a transverse mode in the nth 
lead, while "in" and "out" indicate whether the wave in 
that mode is ingoing or outgoing. As these operators 
obey fermionic commutation relations, one has 



c.-;:: (— e) o-^"' c 



c;in 



i-e) 



(q)T li;in 



(16) 



with a similar relation for outgoing waves. The trans- 
pose in the second term is due to the fact that we com- 
muted the hole operators to ensure normal ordering. We 
then use the scattering matrix to write outgoing oper- 
ators in terms of incoming ones. Contributions com- 



ing from the first term in Eq. ( 16 1 cancel each other 



We find that the operator which gives the spin-current 
along axis a in the electron sector is a a [as in Eq. (13)], 



while it is — cr^ in the hole sector. Recalling that we 



use the convention in Eqs. (2j3), we must also rotate 
the spin-current operator in the hole sector. It becomes 
_^(a)^(a)T^fe) ^ rpj^^g ^j^^ig conven- 

tion, we can write this spin-current operator compactly 
as n'^'^+^a^y) which works for both electrons (/i — 1) 
and holes (/i = —1). From here on, quasiparticle indices 
/i, = +l(e), — l(/i) when they appear as prefactors. 

This calculation in terms of creation and annihilation 
operators for electrons and holes gives us the scatter- 
ing matrix formula that we desire. Assuming that the 
number of transport channels Ni is the same for each 
quasiparticle species, Eqs. (12) is replaced by 



J,. = 



de - 



dl 
de 



J,/3 



(17a) 



E m^^ph - Em- (.)]4« - E M T{r^^"°) (6) 5,p e{T, - To)/To , (17b) 



3-P 



where the integrals now go over a range of positive ex- 
citation energie^and we defined = e{Vj — 14c), i-e. 
voltages are measured from the superconducting voltage 
14c = /^sc/e. We also introduced the spin-dependent. 



quasi-particle resolved transmission coefficients 



^ ^n, r(§A'-)t^(")sMv(«l (18) 
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where is the block of the S-matrix corresponding to 
the transmission of a quasiparticle of type z/ = e, /i in 
lead J to a /i-quasiparticle in lead i. 

The main novelty brought about by superconductivity 
is that the elements of the Onsager matrix now depend 
on Andreev processes via hybrid transmission coefficients 
rj(eh,a0) tj(he,ai}) ^ ^j^jj,]-^ contribute differently to heat 
versus electric and spin currents — see in particular the 
last terms in Eqs. (17a) and (17b). From Eq. ([s]) one 
obtains 



minus signs every time a spin is measured, (iii) exact 
PHS leads to the disappearance of thermoelectric and 
spin caloritronic effects, (iv) at half-filling, exact SLS 
leads to the disappearance of thermoelectric but not spin 
caloritronic effects. 

That thermoelectric and spin caloritronic effects vanish 
in the presence of PHS directly follows from Eq. ( 18 ) 



that transmission coefficients satisfy J^j"^-"^) = 

when PHS is strictly enforced. This gives in particular 



(/ji/;a/3) 



(/ii/:Q0) 



which, together with 



e) = (-l)"°+"''Tj7^^"V3^, e) , (19) Eq. directly gives = T'-^^'i-K) 



_ -p(aO), 



0. 



which extends Eq. \1A\ to include superconductivity. 



Eq. ( 19 ) applies to any hybrid system, regardless of 
whether PHS is present or not. If additionally, the sys- 
tem has unbroken PHS, then the scattering matrix obeys 
Eq. ([t]), i.e. 2>^"'{'K) = ^mu a'^y\$f''' {-K)]* a'^y\ where, as 
before, fiji' — +l{e), —l{h). (Ref. [7] has this formula 
for SRS, where S* commutes with cr^^^). We substitute 
this into Eq. ([l9l, and then substitute cr^^'V^^V^^^ = 



(— 1)"°(T^"'"'". Observing that the trace is invariant un- 
der the transpose of its argument, we find that PHS gives 



7 



ij 



= (-l)"°+"^Tj^^''^"^(-:n,e) 



{-l)^c,+n,jim0<^)^_j^^^y (20) 



ONSAGER RELATIONS 



Eqs. (12), (141, (17), (19l and (20) are all we need 



to derive reciprocity relations between the coefficients of 
the Onsager matrix in Eq. (11). Tables III IV and [V] 
provide a complete list of all Onsager reciprocity rela- 
tions for coupled electric, thermoelectric and spin trans- 
port in single-particle Hamiltonian systems. The On- 
sager relations which can be derived from microreversibil- 
ity are divided into two sets. Firstly, Table [Till gives the 
Peltier/Seebeck relatio ns, between coefficients F^^''^ and 



B^P°l Secondly Table IV gives the reciprocity relations 



for conductances G[J^\ ^i'j^'- example, we note 

that for both the Wigner-Dyson and chiral orthogonal 



classes, the presence of SRS imposes 
while TRS gives T; 



(q/3) _ rr-COO) 



(a/3) 



T' 



(/3a) 



Therefore. 



■'ij "a/3, 

when both 



symmetries are present in those classes, X^"^"^ = Xj^ , 
for X = E, T, B and G. In addition there are those 
Onsager relations which can be derived from either the 
conservation of quasiparticle species (absence of Andreev 
processes turning e into h, and vice- versa), or from the 
presence of PHS or SLS. They are listed in Table [V| 

Some important features are that (i) in multitermi- 
nal devices one needs to consider conductance, Seebeck 
and Peltier matrices, and the reciprocity relations re- 
quire to take their transpose, the latter operation being 
tantamount to momentum inversion as required by mi- 
croreversibility, (ii) spin transport introduces additional 



.(a/3) 



The vanishing of thermoelectric effects with PHS is 
reminiscent of Mott's relation, giving that the Seebeck 
coefficient is proportional to the derivative of the conduc- 
tance at the Fermi energy - the latter vanishes in PHS 
systems. Still, hybrid normal metallic/superconducting 
systems often exhibit larger thermoelectric effects than 
their purely metallic counterpart, which typically hap- 
pens in the crossover regime between Altland-Zirnbauer 
and Wigner-Dyson symmetry classes. For the crossover 
systems described in Table [TTl thermoelectric effects can 
be quite larg^i^HlIl. 

We close this section with two comments on SLS at 
half-filling, when the chemical potential is at zero energy. 
Systems in the chiral symmetry classes have transmis- 
sion coefficients with extra symmetries given in Eq. ( 15 1. 



The latter have important consequences for the symme- 
try of transport, if the trace over the sublattice index 
in Eq. ( 13 ) involves only pairs of sublattice sites, i.e. 
when SLS is not broken by the terminals. When this 
is the case, the firs t an d second equalities in Eq. (151, 
together with Eqs. (flil, give G[f\:K) = G^-f"^(J{) and 



Gffi^) = (-l)"°+"''Gff (-?£) respectively, where we 
recall that no = and n^.y^z — 1- We obtain identical 
results for S.-™'', and thus conclude that 



"(00)/ 



(21) 



We see that charge-conductance, spin-conductances and 
thermal conductance are even in external fields, 5f , irre- 
spective of how many terminals the device has. This is in 
contrast to normal- metallic systems without SLS, where 
only two-terminal devices have conductances even in J{. 
In contrast, spin-to-charge and charge-to-spin conversion 
are strictly odd in 3-C, irrespective of how many terminals 
the device has. 

Turning to thermoelectric and spin caloritronic ef- 



— Fj"^'(Jf), while the second equality in Eq. (151 gives 

us the usual relation b'i^°\:K)To = (-l)"^r)^(- Jf). 
Thus we can conclude that 



fects, the first equahty in Eq. (|l5l) gives B'^f°\j{)To 



S'f )(J{) = Bl^">i~^ for f3 e {x,y,z}, (22a) 



?(/90) 



Sj°)(3{) =0, 



(22b) 
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Symmetry class 


Seebeck-Peltier Onsager relations 
from microreversibility 


Wigner- Dyson 


A (unitary) 
AI (orthog.) 
All (sympl.) 


Bl^To = r'f cc So, 


Chiral 


Alll(unitary) 
RDT ('nrtVinp- 
CII (sympl.) 


Blf'>(J{)To = i-ir^rf/H-K) 

slf^'ro = r'f oc So, 
^r'To = (-i)"^r^f' 


Altland-Zirnbauer 


D 
C 

Dili 
CI 


Bl^'>\ji)To = rf/\-%) 

sr'ro = (-i)"^r(f' 
srro = r(f' ocv 



Table III: The Onsager reciprocity relations arising from microreversibility and which involve the Peltier and spin-Peltier 
matrix elements, r^'^''', and the Seebeck and spin-Seebeck matrix elements, B'^^'^K These relations, combined with those due 
to PHS or SLS in Table [y] give the complete set of Onsager relations for each symmetry class. 



Symmetry class 


Onsager relations between conductances, A = G, H 
from microreversibility 


Wigner-Dyson 


A (unitary) 
AI (orthog.) 
All (sympl.) 


X[f\-K) = (-l)"°+"'3x]f">(-M) 

-A -J - oc da, 


Chiral 


A Ill(unitary) 
BDI (orthog.) 
CII (sympl.) 


X^f'^M) = (-l)"°+"^X]f"'(-5{) 
Xij - X^- OC da. 


Altland-Zirnbauer 


D 

C 

Dili 

CI 


xlf\M) = (-l)"°+"3X]f">(-M) 

A(;«(J{)=A]f"'(-Jf) 
Xij - OC da. 



Table IV: The Onsager reciprocity relations arising from microreversibility and which involve the electrical and spin-dependent 
conductances G^"^^ and the heat conductance S''"'', These relations, combined with those due to PHS or SLS in Table |v] give 
the complete set of Onsager relations for each symmetry class. 



with identical relations for F 



(0,3) 



Additionally, any sys- 



tem with PHS has no thermoelectric nor spin caloritronic 
response. Looking at Table |T] we see that only the AIII 
symmetry class has SLS without PHS. Thus in this sym- 
metry class, the spin-Seebeck and spin-Peltier coefficients 
are even functions of the external field, IK, while the usual 
Seebeck and Peltier coefficients vanish identically. 



We stress, however, that the analysis leading to 



Eqs. (21) and (22) holds only at half-filling, when the 



Fermi function in Eq. ( 12 ) is symmetric around e = 0, 



and, perhaps physically more important, when the ter- 
minals do not break SLS. This requires leads to be con- 
nected with equal strength to both sublattice sites in each 
unit cell. 



VI. EXAMPLES OF RECIPROCITY 
RELATIONS IN SPINTRONICS AND SPIN 
CALORITRONICS 

A. Spin Hall and inverse spin Hall effects 

As a first example o f the r eciprocities we derived, 
we discuss th e spin HalP^'^^'^^^'Sfil and the inverse spin 
Hall22126H30l53] effects. The two effects are sketched in 
Fig. |2] In the spin Hall effect, Fig. [2^, one passes an 
electric current between terminals 1 and 2 and measures 
the spin current between terminals 3 and 4. The voltages 
at terminals 3 and 4 are set such that no current flows 
through them on time average. In the limit of large and 
identical number of channels in each terminal, iV ^ 1, 
the voltages V3 and V4 lie almost exactly in the middle 
between Vi and V2, VzA — (^1 + ^2)/2 for balfistic sys- 
temPl. We assume that this is the case here, and set 
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Symmetry class 


0|J(:3d<3.1 dLlUiLlUIicli I t^lclLlUiiS 


Wigner 
— Dy^son 


A (unitary) 

A T (' nrthnp' ^ 

All (sympl.) 


n(00) _ n(00) 

ij ji 
^(00) ^ ^(00) 


Chiral 

(half-filling only) 


Alll(unitary) 
BDI (orthog.) 


(M) = (-M) for X = {G, H}, S'f (M) = B^f °' (-J{), Bff (J{) = 

d(/30) _ p(0/3) _ 
n(/30) _ p(0/3) _ „ 


Altland 
-Zirnbauer 


D 
C 

Dili 
CI 


-Bij (^) = Mj (W) = 

B(f)(jt) = r(°«(w) = o 

d(/30) _ p(0/3) _ r, 

o(/30) _ p(0/3) _ „ 
ij ij 



Table V: Additional reciprocity relations induced by the conservation of each quasiparticle species (absence of Andreev reflection 
from e to h), or by the presence of PHS or SLS, to be added to the Onsager relations in Tables III and IV The additional 
relation between B and its transpose (and an identical one, not listed here, between F and its transpose) in the Wigner-Dyson 
classes was first noticed in Ref. [44], which allows one to express the Seebeck/Peltier relation of Table III in two different but 
equivalent ways. Thermoelectric and spin caloritronic effects disappear identically in the presence of PHS. The relations in the 
chiral classes correspond to transport at half-filling, E — 0. Outside of this regime, chiral systems have the same relations as 
in the (corresponding) Wigner-Dyson classes. 



Vi = V/2, V2 = -V/2, V3,4 = 0. 

The presence of spin-orbit coupling inside the system 
generates a spin current flowing through the transverse 
terminals. When all terminals are at zero temperature, 
these currents are given by 



Together with Eq. (14), Eq. (25a) gives GsHc — GisHc- 
The reciprocity between direct and inverse spin Hall con- 
ductances is exact and does not require sample averaging, 
as sometimes claimec 



2h 
2h 



\-'31 ■'32 )j 



(aO)x 



■'42 ) ■ 



(23a) 
(23b) 



In the inverse spin Hall effect. Fig. [2|d, there is no volt- 
age bias, but instead terminals 3 and 4 have opposite 
spin accumulations. In an idealized situation they will 
be ±/i^"^/2. Spin-orbit coupling converts this spin accu- 
mulation into a transverse electric current. The currents 
in terminals 1 and 2 read 



r(o) 



r(o) 



2h 



-{7\ 



13 



7\ 



' 2h ^ ■ 



23 



<Oa)x 
14 ) ■ 

^-(0a)^ 
■'24 I ■ 



(24a) 
(24b) 



In both cases, the Hall part of the currents, flowing be- 
tween 3 and 4 in the case of the spin Hall effect and 
betwee 1 and 2 in the case of the inverse spin Hall ef- 
fect, is given by the difference in the two currents. We 
define spin Hall and inverse spin Hall conductances as 



;-(") 
^sHe 

GisHe 



r(«) 
^3 

»/e. 



GsHeV and I, 



(0) 

isHe 



(0) 



One obtains 



GsHo — 
GisHe — 



2/1^ 31 



•^32 ■'41 



7 



e 
'2}i 



(■''23 



(Oa) 



7: 



(Oq) 
24 



7- 



iOa) 
13 



<0q) 
14 



r(0) _ 
^1 — 



(25a) 
(25b) 



Reciprocity between spin injection and 
magnetoelectric spin currents 



For a spin index /3 = 0, Eqs. (14) and (19) establish 



the reciprocity between magnetoeletric effects generat- 
ing spin currents from electric voltage biases and spin 
injection from spin accumulations in the terminals, a 
special case of which is the above-discussed spin Hall ef- 
fect /inverse spin Hall effect reciprocity. In the presence of 
TRS, it has already been observed that one consequence 
of Eq. ([5]) is that no spin current can be magnetoelectri- 
cally generated in a two-terminal device if the exit lead 
carries a single (spin-degenerate) transport channel. The 
reciprocity relations of Eqs. ( 14 ) and ( 19 ) further impose 



that a spin injection from such a terminal is incapable 
of generating an electric current, unless one goes to the 
nonlinear regimd^. This seems not to have been noted 
so far. 



C. Spin Seebeck and spin Peltier coefficients in 
two-terminal geometries 

In two-terminal geometries, the electric conductance 
is symmetric in TRS breaking fields, which follows 
from current conservation or gauge invariance, together 
with the symmetry of electric reflection coefficients. 



G'^,\'K) = G^Pi-'K) (see e.g. Ref. gT]). Including 
spin-transport, the unitarity of the scattering matrix fur- 
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Figure 2: Sketch of a four-terminal spin-Hall (a) and inverse spin Hall (b) experiment, (a) In the spin Hall effect, an electric 
current (blue arrow) generates transverse spin currents (red) via the action of spin-orbit coupling, (b) In the inverse spin Hall 
effect spin accumulations inject a spin current (blue arrow) which, in the presence of spin-orbit coupling, generates a transverse 
electric current (red). 



ther results in spin-current conservation and generalized 
gauge invariance, 



J2 {msoo^s^j - t; 



t: 



(Oa) 
ij 

(qO) 



0, 



= 0, 



(26) 



under the assumption that the number of transport chan- 
nels coupling the system to external reservoirs is spin- 
independent. In two-terminal geometries this gives 



B 



m 
11 



B 



m 

12 



B. 



m 

21 



r(o/J) 
^ 11 



^21 — ^ 12 



r 



(PO) 
22 

(0/3) _ 



0, 



0. 



(27) 



However, unlike for the charge conductance, the thermo- 
electric reflection coefficients can have both a symmetric 
and an antisymmetric component. This is directly seen 
from the expression 



de 



dl 
de 



-to 



(28) 



for the spin Seebeck reflection coefficient. For example 
for (3 = z, the spin-dependent transmission coefficient in 
the integrand reads 



•' ii 
•' s 

7(20) 



s •'a ' 



T, 



T, 



(29a) 
(29b) 
(29c) 



which, from Eq. i5 



7e^^\—3-C), and antisymmetric, 'ji^°'(3{) 



has both symmetric, 7i^^\3-C) — 



components. 

An interesting example is provided by a two-terminal 
system with a well-defined spin quantization axis. This 
is the case, for example, for a system without spin-orbit 
coupling in a uniform Zeeman field, for two-dimensional 
systems with both Rashba and Dresselhaus spin-orbit 



interactions of equal strength^SZI^ or for a system with 
pure / • s spin-orbit coupling. Without loss of gener- 
ality we define the spin quantization axis as the ^-axis. 
Then § commutes with a(^), i.e. it is diagonal in spin 
space. From Eq. (fisl), we find that T[f ^(J{) = T^f ^(:H;) 



and t|^'(M) = when a = x,y. Combining this with 
Eqs. Q, we have T^f (J{) = 'j'^2'\o<) = T^f (^K) ^ 
1^21 ^(Di). Thus 



B[f\:K)To = TY2>{:K), 



_ p(Oz), 



(30) 



withr^2"''(^) = Owhen a = x,y. Next we recall that the 
Seebeck-Peltier Onsager relations contain an extra minus 
sign for spin caloritronic effects compared to usual ther- 
moelectric effects (see Table III). This extra minus sign 



means that 5^2°'' ^'^'^ ^12' ^-re odd in IK, while B'"y^' and 
r^2'''' E^re even in 'K. Thus any two-terminal system with 
a spin-quantization axis will have spin-Seebeck and spin- 
Peltier effects which are odd functions of TRS breaking 
fields, while the normal Seebeck and Peltier effects are 
even function of those fields. 



(00) 



VII. EXAMPLES OF RECIPROCITY 
RELATIONS IN THERMOELECTRICITY WITH 
HYBRID SYSTEMS 

Thermoelectric effects in the presence of superconduc- 
tivity, in particular the thermopower S = — _b(oo)/G'(°'') 
and thermal conductance s'""^ have attracte d quite 
some experimentaP^lti^l and theoretical interesiP *^^ l ^^ l ^^ l 
However, the exact form that the Seebeck-Peltier On- 
sager reciprocity relation takes has never been clari- 
fied, despite the fact that two-terminal devices with su- 
perconductors usually exhibit odd Seebeck coefficients 
S{'K) = -S{-'K), in stark contrast with Mott's re- 
latiorPSi. Mott's relation between the thermopower of 
metallic systems at low temperature and the energy 



11 



superconducting loop 



Lead 1 
(normal) 




Lead 2 
normal) 



Figure 3: An Andreev interferometer: two-terminal hybrid 
system with two normal metallic/superconducting contacts 
with tunable superconducting phase difference (j}sc- The lat- 
ter, together with a systematic delay (indicated by the extra 
length £) towards one of the superconducting contacts can 
break PHS and generate finite thermoelectric effect^^. The 
particular symmetry class in which the system falls is given by 
how much magnetic field 3f and how much spin-orbit coupling 
there is in the dot, see Table |l] 



derivative of the conductance at the Fermi energy, reads 



5* = -- 



3e 



(31) 



and thereby indicates that S should be even in This 
evenness of S is confirmed by the scattering theory for 
metallic systemsA^. In this section we provide examples 
clarifying this issue using scattering theory to show that 
S can have any symmetry under !K — > — Jf when super- 
conductors are present. 



A. Seebeck— Peltier reciprocity relation 

Andreev scattering strongly influences the Seebeck- 
Peltier reciprocity relation between F and B coefficients. 
Comparison of the last terms in Eqs. (17a I and (17b) 
shows that the ee and hh terms in F and B have the 
same sign, while the eh and he terms acquire a relative 
minus sign. This breaks one of the Onsager relations 
between Peltier and Seebeck coefficients. For metallic 



systems, one has both i?-™^(IK)To 



b^T\k)To 



= F 



(00) 



and 

however, with superconduc- 



.(00) 



ity, J2j V '^iji^^ ^) = ^'^'^ assuming that the number 
= and both F- and B-coefficients vanish of transport channels depends neither on the quasi- 

particle type nor on the magnetic field, we readily obtain 



'ity, only 4™)(J{)To = F^f holds. 
When PHS strictly holds, however, J2u vT'i'^^'"^'' 

identically, regardless of the temperature. However, in- 
teresting thermoelectric effects appear in hybrid systems 
when PHS is broken. Focusing on a two-terminal geom- 
etry, as depicted in Fig. |3j Eq. (17) can be rewritten in 
the form 



(32) 



which depends only on the voltage and temperature dif- 
ferences between the two normal reservoirs. The two- 



terminal thermoelectric coefficients are given by 



G = 



^(00) 



^(00) 

"11 



^(00) 



(00) 



i(00)\ 



^(00) 

'^ll 



G 



(00) 
22 



^(00) 



B 



B 



(00) 

11 



^(00) 
^(00) 



g: 



(00) 



Gi2°' 

(OO)wp(OO) 



22 



{G^^>^>^+G^^!^'>){B[l'>^+B^f^) 



+ "-^21 
?(00)^, 



.(00) 
11 



^(00) 

^r(oo) _ 



G, 



(00) 



^(00) 



Fir)(Gir+GD 



^(00) 
(00)\ 



'-'11 + '-'2 



r^(oo) ^(00) 
'-'12 +'-'21 



(33a) 



(33b) 



(33c) 



(33d) 



in terms of the coefficients X'™' (X — G, B, F, ^) defined 



by Eqs. (11) and (17 1 



It is then straightforward to see that the reciprocity 
relations read specifically 



g(:k) 

5(Jf) To 



G(-J{), 
S(-5{), 

f(-:k) . 



(34a) 
(34b) 
(34c) 



In particular the presence of superconductivity forces one 
to inve rt the sign of the TRS breaking field in the relation 
of Eq. (34c) between Seebeck and Peltier coefficients. 



B. Symmetry of the thermopower 

The symmetry of the two-terminal thermopower, S — 
__B(oo)/(5(oo) jg jjQi- specified in the presence of super- 
conduct ivitjff^. The Seebeck coefficients read 



2e 

Wo 



dsi~dj)e 



""^(£,J{)-T(;''^°°)(£,5£)].(35) 



-t: 



(/ie;00) , 



From this expression we see that thermoelectric effects 



vanish, B 



(00) 



0, if PHS is enforced; we thus con- 



sider this equation in the absence of PHS. From Eq. ([s]) 



we know that T, 



(w;00) 



T 



■(eft;00) 



(e,J{) 



_ rj-(/J/j;00) j-^ 



M), while 
Together with unitar- 



that bJ°^(:K) = -Bin (IK) + 
even and an odd component. 



(00), 




2e 



x[T5^^°")(e,Jf)-Tg^'^^°°)(£,J£)] 
2e 



Ae{-dJ)e 



X 



rrp(e;i;00) 



(36a) 



(36b) 
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where Bi°°li-:K) = Bi°S(Jf) and = 
— i?^™^(5f). In the absence of Andreev scattering, 

S(™)(3-C) = i3even(3^) is strictly even in two-terminal ge- 
ometries, however Andreev scattering gives rise to an odd 
component. The asymmetric Andreev interferometers 
considered in Ref. [T^ were devised to render B^^^ {3i) fi- 
nite on mesoscopic average, which led to an antisymmet- 
ric thermopower in such systems. There are currently 
no known hybrid systems which have a finite-average 
-Seven (3^). Recent theoretical works pointed out asym- 
metries in the thermopower of metallic systems in the 
presence of inelastic scattering, which is of interest be- 
cause asymmetric thermopower may lead to more efh- 
cient thermal engineJ^^^. Hybrid systems are examples 
of systems with purely elastic scattering and antisym- 
metric thermopower. 



Onset of thermoelectric effects upon breaking of 
PHS 



Thermoelectric effects vanish identically in all Altland- 
Zirnbauer symmetry classes because of PHS. However in 
physical systems PHS is often at least partially broken, 
leading to finite thermoelectric effects. Here we show 
that the symmetry of such thermoelectric effects is subtly 
dependent on how PHS symmetry is broken. 

To that end we consider the Andreev interferometer 
shown in Fig. |3] A two-terminal chaotic ballistic or dis- 
ordered diffusive quantum dot is connected to a super- 
conducting loop via two contacts. The superconducting 
phase difference at the two contacts can be tuned by a 
magnetic flux piercing the loop. There are two impor- 
tant time scales in the system, (i) the typical time TAndr 
between two consecutive Andreev reflections at the su- 
perconducting contact, and (ii) the escape time Tcse to 
one of the normal leads. We additionally choose a special 
geometry where the average time to reach one of the two 
superconducting contacts from one of the normal leads 
is longer - this is achieved by an extra ballistic "neck" 
of length £ between the cavity and the superconducting 
contact (see Fig. [3). Because of the neck, quasiparticles 
need an additional time delay St = £/vf to reach the 
left superconducting contact from a normal lead. To- 
gether with this time delay, a magnetic flux piercing the 
superconducting loop and making the superconducting 
phase difference (f>sc finite al so bre aks PHS, thereby turn- 
ing thermoelectric effects oi i^QI^^I . 

Formally, PHS requires that TAndr — > 0, which practi- 
cally means that TAndr has to be smaller than any other 
time scale and any other inverse energy scale. When this 
is not the case, transport processes without any Andreev 
reflection exist, giving contributions to the conductance 
that fluctuate randomly in energy around the Fermi en- 
ergy. This breaks PHS and leads for in stance to finite, 
albeit relatively weak, thermopowei'S21I2! More generally, 
breaking PHS can be achieved in three different ways. 



(i) rendering escape into the normal leads faster (for 
instance by widening the normal leads), until 



TAndr i 



(ii) raising the temperature until (h-BT) ^ ^ TAndr, or 

(iii) changing the flux through the superconducting loop 
so that 7^ Oi ""j when the neck length £ is finite. 



In case (i), a significant proportion of quasiparticles go 
from one normal lead to another without Andreev reflec- 
tion. Then contributions to •J^^j'^^"'^^ (3^) which arise from 
processes without Andreev reflections will start to dom- 
inate thermoelectric transport, meaning Seven 3> B'^^J^ 
[as defined in Eq. (36)]. Thus, thermoelectric effects ac- 



quire the same symmetry as systems without SC con- 
tacts, i.e. they become predominantly even. 

The situation is more complicated in case (ii), where 
both t|^'"''"^^(J{) and 7'if'-"'^\^) have similar magni- 
tude. In the absence of a neck, ^ = 0, thermoelectric ef- 
fects vanis h on a verage and are dominated by mesoscopic 
fluct uat ion^i^liiJ. An analysis of these mesoscopic fluctu- 
ations analogous to that in Ref. |T0] shows that there 
is no correlation between ''(?£) and Bj^^\—^), so 
that the thermoelectric effects have no particular sym- 
metry _beyond the generic Onsager reciprocities given in 
Table 



III 

TOO) 



(00) 
even 



In particular, for a two terminal device B, 
and B^^J are independent random variables with the 
same variance. Thus for a given Andreev interferome- 
ter (given disorder or cavity shape) either quantity could 
be positive or negative, and either could have a larger 
magnitude than the other. 

Finally in case (iii), the physics changes completely. 
Due to the presence of a finite-sized neck, £ 0, and 
superconducting phase difference (jjgc 7^ 0,7r, the system 
develops a large average thermopower which is an odd 
function of the flux (/(gj^'^, with a much smaller even 
component coming from mesoscopic fluctuationJ^Sl 

In summary depending on how particle-hole symme- 
try is broken, one gets a thermopower which is predomi- 
nantly even in 5{ [case (i)], predominantly odd in Jf [case 
(iii)], or which has no particular symmetry [case (ii)]. 



VIII. CONCLUSIONS 

We have derived a complete list of reciprocity rela- 
tions for coupled electric, spin, thermoelectric and spin 
caloritronic transport effects in all ten symmetry classes 
for single-particle Hamiltonian systems. Several of these 
relations appeared in one way or another in earlier works, 
and the main novelties we found are (i) reciprocities in 
spintronics and spin caloritronics pick a number of ad- 
ditional minus signs reflecting spin current injection and 
measurement, (ii) a number of special relations have been 
listed in Table |Vj which exist only in specific symmetry 
classes, (iii) we clarified the exact form of Onsager rela- 
tions in the presence of superconductivity, and (iv) we 
derived all Onsager relations for transport in spintronics 
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ALL SYSTEMS 



'{n)To = {-irT'"'>{-n) 



No Andreev Refl. 



Time-reverse sym. TRS= 
Magnetic-field H = 



±1. 



Spin— rotation sym. TRS= +1. 



al3 



Bir - 



Figure 4: Venn diagram summarizing the Onsager relations we derived. Here a, /3 £ {0, x, y, z}, no = and = riy = Uz — 1. 
Any system obeys the relations given in all boxes it is within. For instance the Onsager relations for a generic system with time- 
reversal symmetry (TRS) read G^J*^' = (^_i'^n^+n„ ^ipo,) 3^^°^ To = (-l)"«rjf \ since TRS requires that any external 
magnetic field J{ = 0. If that system also had SRS, then the only nonzero coefficients are those with repeated upper-indices 

If the system contains no Andreev reflection (i.e. 



(/3/3) 



(a = /3 for G and /3 = for B) for which G^^ 
no superconductors) , one additionally has'** B, 



derived in the text are summarized in Table IVl 



G'f) and Sf^To 



-.(00) 



-^I'i"' (^-^)- case of PHS and SLS, the extra symmetry relations 



and spin caloritronics in the presence of superconduc- 
tivity. We present a pictorial summary of the Onsager 
reciprocity relations we derived in Fig. [4] 

Generally speaking, our investigations of the specific 
reciprocities shown in Table |V] allowed us to clarify 
the form that the Seebeck-Peltier relations take in the 
presence of superconductivity. While the two relations, 



B, 



(00), 



= (:K) and Bl^'^>{K)To - 
exist in purely metallic systems, only one of these two 
Onsager relations survives in the presence of supercon- 
ductivity, that being B^^°\'K)To = r5™^(-J{). 



(00), 
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